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Abstract 

It is shown that quantum vacuum fluctuations give rise to a cur- 
vature of space-time equivalent to a cosmological constant, that is a 
homogeneous energy density p and pressure p fulfilling -p = p > 0. 
The fact that the fluctuations produce curvature, even if the vacuum 
expectation of the energy vanishes, is a consequence of the non-linear 
character of the Einstein equation. A calculation is made, involving 
plausible hypotheses within quantized gravity, which establishes a re- 
lation between the two-point correlation of the vacuum fluctuations 
and the space-time curvature. Arguments are given which suggest 
that the density p might be of order the "dark energy" density cur- 
rently assumed to explain the observed accelerated expansion of the 
universe. 

Keywords Dark energy. Vacuum fluctuations. Quantum gravity 
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1 Introduction 

Recent astronomical observations, in particular the study of type la super- 
novae, anisotropies in the cosmic background radiation and matter power 
spectra inferred from large galaxy surveys, have improved our knowledge of 
the universe giving rise to a precision cosmology. The new data are compat- 
ible with the universe having a Friedmann - Robertson - Walker metric with 
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flat spatial slices pQ of the form 



ds^ = -dt^ + a{tf {dr^ + r^dn) . 



(1) 



The time-dependent parameter a{t) is related, at present time Iq, to the 
measurable Hubble constant, Hq, and decceleration parameter, go? via 
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The observations also provide information about the evolution of the param- 
eter a{t) and the distribution of matter in the past. 

The available knowledge is summarized in the ACDM model. In it bary- 
onic matter density, p^, represents about 4.6% of the matter content of the 
universe while two hypothetical ingredients named cold dark matter (CDM) 
and dark energy (DE) contribute with densities p^j^j ~ 23%, and p^^ ~ 
73% respectively [2]. The said densities are related to the metric eq.(IT]) via 
the Friedmann equations, derived from general relativity, giving the following 
relations [3] 
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where I neglect small effects of radiation and matter pressure. The baryonic 
density p^ is well known from the measured abundances of light chemical 
elements, which allows calculating pj^^ and p^^^ from the empirical quantities 
Hq and go via comparison of eqs.(I3]) and The values obtained by this 
method agree with data from other observations. For instance cold dark 
matter, in an amount compatible with Pdm, is needed in order to explain 
the observed (almost flat) rotation curves in galaxies. However the nature of 
dark matter and dark energy remain open problems. 

Usually dark matter is assumed to derive from exotic particles, not yet 
discovered, and dark energy is identified with a cosmological constant, A. In 
any case the ACDM model rest upon the assumption that general relativity 
(GR) is indeed the correct theory of gravity. However it is conceivable that 
both cosmic speed up and dark matter represent signals of a breakdown of 
GR. For instance we might consider the possibility that the Hilbert - Einstein 
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Lagrangian, linear in the Ricci scalar R, should be generalized to become a 
function f{R). This is the underlying philosophy of what is referred to as 
f{R) gravity theoryjl]. Indeed the cosmological constant corresponds to a 
particular choice of f{R) where a constant A is added to the Ricci scalar 
R, although this does not give any hint about the value of pj^^. The theory 
of f{R) gravity provides sufficient freedom to accomodate also dark matter. 
For instance it allows good fits to the rotation curves in galaxies[S], which 
therefore might be explained as a curvature effect. Actually f{R) gravity 
may be further generalized by including other scalars, like R^^R'^'^, in the 
Hilbert - Einstein Lagrangian [6]. In this paper I shall not deal with dark 
matter, but only with dark energy, so that I will take the density pj^j^j as an 
empirical datum, without any discussion about its possible nature. 

Many proposals have been made for the origin of dark energy (for a 
review see Copeland^). As said above the most popular one is to identify 
it with a cosmological constant or, what is equivalent in practice, to assume 
that it derives from the quantum vacuum. Indeed the term pj-)^ in eqs.(l3]) 
might be interpreted as coming from a vacuum energy whose pressure fulfils 
Pde = ~Pde^ equality appropriate for the vacuum (in Minkowski space, 
or when the space-time curvature is small) because it is invariant under 
Lorentz transformations. A problem appears however when one attempts 
to estimate the value of pj^^. For instance if the dark energy is due to the 
interplay between quantum mechanics and gravity, it may seem that it should 
be either strictly zero or of order Planck's density, that is 

Pde--§^^^ 10'^ kg/ml (4) 
In sharp contrast the observations lead to 

Pj,e ^ 10-'^ kg/m^ (5) 

a value departing from eq.(|l]) by about 123 orders of magnitude. This strong 
disagreement gives rise to the "cosmological constant problem" |8] . 

In this paper I will explain why the value of pj^^ can be much smaller 
than eq.(jl]) even if the cosmological constant really derives from the quan- 
tum vacuum and I shall do that without departing from standard general 
relativity. In fact I shall prove, modulo a few plausible assumptions within 
quantum gravity, that vacuum fluctuations give rise to a curvature of space- 
time fully equivalent to the one produced by a cosmological constant, even 
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if the vacuum expectation of the density of quantum fields is strictly zero. 
The argument is as follows. 

For a small enough region of the universe around us, but large in com- 
parison with typical distances between galaxies, the space-time metric given 
by eqs.(IT]) may be rewritten, near present time, using new coordinates as 
follows [3] 
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(6) 



where terms of order O (r*^) have been neglected and it is assumed that the 
(slow) change of the metric coefficients with time may be ignored. This met- 
ric is Minkowskian for small r, which makes the calculations more simple 
than using eq.([l]). In this paper I will calculate the coefficients grr and gu of 
eq.(l6]) as coming from the combined action of cold matter, having homoge- 
neous density + pj^jyj (at the large scale), plus the effect of the vacuum 
fluctuations. In my approach the Fnedmann eqs.l^ are not valid because 
they were derived under the assumption that the space-time curvature, re- 
sulting in the metric eq.Q , comes from a mixture of three fluids with total 
density P^ + Pdm + Pde ^^^1 total pressure —pjj^- Here I will assume only two 
fluid with total density p^ + pj^j^j and neghgible pressure. The calculation 
leads to the following relations 

^ {Pb {t) + Pdm {t)) + ^fiuct r'^ + O (r^) , 
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where A jiuct is a constant parameter with dimentions of inverse length squared. 
It is explicitly calculated from the two-point correlation of the vacuum fluc- 
tuations (see below). The net result is that the fluctuations produce the same 
effect on the space-time curvature as a cosmological constant. Thus eq.([7]) 
may be written in a form similar to Friedman's eqs.(]3]) provided we deflne a 
new quantity, p^i^;, with dimensions of density as follows 
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But pj^^ is not any actual density, but a parameter taking into account the ef- 
fect of the quantum vacuum fluctuations on space-time curvature. My deriva- 
tion does not provide a precise value of but it strongly suggests that it 
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is far smaller than eq. . The difficulty for getting the value is that the 
two-point correlation of the vacuum fluctuations is not known. If future 
calculations along this line provide a value of pjj^ in agreement with eq.® 
then the universe speed up would be fully explained as due to the quantum 
vacuum fluctuations. 



2 The two-point correlation function of vac- 
uum fluctuations 

The starting point of the work is an idea of Zeldovich[n], who proposed the 
relations 



PoeC -G— = — y.-„X^- (9) 



where m is a typical mass of elementary particles, eq.® being obtained if 
the mass m is 

m ~ 7.6 X 10"^^ kg, 

that is about one third the pion mass. In my opinion the Zeldovich 's relation 
between cosmology and particle physics suggests that dark energy density 
eq.® does not correspond to the mean vacuum energy, which is likely zero, 
but to (small) gravity effects associated to the quantum vacuum fluctuations. 

The second eq.® suggests that PdeC^ may have the magnitude of the 
gravitational energy of quantum vacuum fluctuations. This may be seen more 
explicitly using a semiclassical Newtonian theory of gravity, that is taking 
matter as quantized but the gravitational field as classical. Thus we may 
calculate the vacuum expectation of the gravitational energy associated to a 
sphere of radius R, that is 
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(10) 

This calculation rests upon the assumption that it is possible to define an 
energy density operator, p (r,t) , of the quantum fields and that the vac- 
uum expectation of that energy is zero but the expectation of the square is 

a2 

not zero. Actually if the expectation of p was also zero there would be no 
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quantum fluctuations at all. That is we must assume 
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(11) 



This being the case, by continuity we expect a two-point correlation func- 
tion which may depend only on the distance |r2 — ri| for equal times in the 
non-relativistic approach leading to eq.f lT(I]) . In a relativistic theory the cor- 
relation should depend on the interval, s, this being the only invariant in 
Minkowski space. The operators p (ri,ti) and p (r2,t2) may not commute 
and it is plausible to define the two-point correlation with the operators 
in symmetrical order, which guarantees that the correlation, C{s), is real, 
without an imaginary part. That is 



C(s) 



vac 



p (ri, ti) p (r2, t2) + p (r2, tg) p {ri,ti] 
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with units c = 1 which I shall use from now on. With this definition the 
correlation might depend on whether the interval s is space-like or time-like, 
in the latter case s being imaginary. 

Taking eq.( IT2|) into account and considering equal times, ti = ^2, eq.( !TO|) 
leads to 



E 



'\r2\<R J\r2\<R 



. C'(|ri-r2|) 



|ri - r2 



(13) 



-G I rf^ri 

|r2|<-R ^0 



47rr — —dr 



-4t:GV I C{r)rdr, 
'o 



where r stands for |ri — r2| and the integral in r has been extended to oo 
because we assume that the radius R is much larger than the range of the 
correlation function C (r) . The result shows that, in Newtonian gravity, the 
existence of fluctuations necessarily implies a gravitational energy associated 
to them, with density 
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The gravitational energy appears in spite of the fact that the vacuum ex- 
pectation of the matter density operator vanishes everywhere, see eq. (fTT|) . 



6 



The fact that semiclassical Newtonian theory predicts that vacuum fluctua- 
tions give rise to a gravitational energy, even if the mean energy density of 
the vacuum is strictly zero, shows that in semiclassical general relativity the 
vacuum fluctuations will produce space-time curvature. Indeed in general 
relativity the concepts of gravitational force and gravitational energy lose 
their meaning and the relevant concept is the curvature of space-time. 

The aim of this paper is to present a calculation using general relativity. 
Thus the goal of our calculation will be to flnd the curvature of space-time in- 
duced by the quantum vacuum fluctuations. The calculation develops an idea 
put forward elsewhere pT]. Before presenting the calculation I shall discuss 
further the subject of the two-point correlations of vacuum energy densities. 
This is necessary because quantum vacuum fluctuations might be seen as ar- 
tifacts of the quantum formalism, without real physical implications, in view 
that in most cases they may be eliminated by using normal ordering of the 
creation and annihilation operators. 

The two-point function C (s) might be calculated in flat (Minkowski) 
space-time from the properties of quantum flelds in vacuum, but making 
a calculation which involves all known flelds would be a formidable task. 
Nevertheless the calculation is straightforward in principle, as shown by the 
derivation which follows of the contribution due to the free electromagnetic 
fleld, which I will perform for illustrative purposes. In quantum theory the 
vacuum expectation of the energy of any unexcited fleld is assumed to be 
zero and this assumption is stated formally by using the normal ordering of 
the operators. For instance in the electromagnetic fleld we have 

{vac\pir,t)\vac) = 0,pir,t)= : ^ ' \ : , (15) 

where normal ordering implies that the energy density operator p contains 
only products of creation, a"*", and annihilation, a, operators of photons of the 
type oa, a^a^ y a^a, all of which give a nil vacuum expectation. Therefore 
the vacuum expectation of the energy density is zero as assumed. However 
the two-point correlation is not zero because the operator 

_ ^ (16) 
contains terms of the type aaa~^a'^ whose vacuum expectation is flnite. 

The two-point correlation function of the free electromagnetic fleld will 

be just the vacuum expectation of eq.( IT6|) . The calculation is straightforward 



using the plane-waves expansions 
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with standard notation {V is a normalization volume and u = |k|). The 
details may be seen in the Appendix and the result is 



The two-point correlation depends on distance and time interval only via 
— as it should. Furthermore it depends on — , so making no 
distinction between time-like and space-like intervals. 

The correlation C (s) , eq.f lTSi) , decreases rapidly at large — but has 
a strong divergence when — )■ t^. It is plausible to assume that quantum 
fields other than the electromagnetic one will give rise to counterterms which 
eliminate the divergence. In particular the sums in k which appear in eqs.( !T7|) 
have been extended to very large values of |k| , but this is physically absurd. 
Indeed the plane waves expansions eqs. f ll7p correspond to assuming an energy 
^fiu per normal mode of the radiation but, for values of fiu larger than 
two electron masses, the electromagnetic field may produce electron-positron 
pairs and the study of the radiation field alone makes no sense. This may be 
also stated saying that at very high frequencies the electromagnetic vacuum is 
polarized, so that the denominators Sir in eq. f[T^ should be replaced by larger 
quantities. For a rigorous treatment we should study the electromagnetic 
field interacting with the electron-positron field, but we should also include all 
other charged particles, and all particles interacting with these via strong of 
weak forces. In summary, a consistent calculation of the correlation function 
should involve all quantum fields (excluding gravity.) In the absence of that 
calculation 1 shall assume that the tail of the correlation function C [s) is 
given by eq.f fTSj) . the photon being the only known massless particle which 
may exist freely, but that at small values of s the function C (s) remains 
finite. If this is the case, it implies that the contributions of particles with 
different masses are not additive, but may cancel each other to some extent. 




(18) 
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Another mechanism able to remove the strong divergence of eq.f llSp will be 
discussed at the end of this section. A simple form to take into account these 
possibilities is to introduce a cut-off, substituting the following for eq.f fTSj) 

where m is an unknown mass. Putting this into eq. ffT^ we get 

4 Gm^ 
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grav 



37r3 fi4 



close to Zeldovich 's proposal eq. ([9]) . It is not easy to estimate the value of the 
cut-off mass m without a detailed knowledge of the correlation function C(s), 
eq.f lT^ . but we may guess that m is of order the masses of the fundamental 
particles of the "standard model" rather than of order Planck's mass. Thus 
a model resting upon Zeldovich 's idea predicts a density pj^^ much closer, 
if not identical, to eq.(l5]) than to eq.(|l]) . 

The dependence of the two-point correlation on the interval s gives rise 
to a paradox which reflects the counterintuitive feature that the correlation 
does not decrease with distance. It may be seen as a straightforward predic- 
tion of relativistic theory, but I shall show in the following that it is really 
counterintuitive. The paradox may be stated as follows. Let us assume, for 
the sake of simplicity, that the density vacuum fluctuation in a point pos- 
sesses discrete values, say with probability Pi, p2 with probability P2, etc. 
Then the mean square fluctuation will be (compare with eq.f lT^ ) 

C(0) = 5^P,pf. (19) 

3 

Now we consider two different points, the correlation of fluctuations being 

^(^) = EEwiPfc, (20) 

3 k 

where pjk is the probability that the vacuum fluctuation in the first point is 
Pj and the fluctuation in the second point it is p^. If the two said points are 
separated by a light-like interval, then s = so that eqs.f fT9|) and fl20l) lead 
to 

j k j j k 
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where the latter equahty follows from well known properties of the probabil- 
ities. Hence it is trivial to derive the equality 

J2J2pi'' + Pk- '^PjPk) = J^J^P^i" (Pj - Pkf = 0- 

j k j k 

For (positive) probabilities this equality can be true only if pjk = for any 
j 7^ k. This means that the fluctuations are strictly correlated in the whole 
light cone of every point. Furthermore, for two arbitrary points, 5*1 and 5*2, 
it is always possible to find another point 5* which is light-like separated 
from each one. In fact, all points in the intersection of the light cones of 
Si and 5*2 do the job. As a consequence for all pairs of points the proba- 
bilities pjk are zero for any j ^ k. which implies that vacuum fluctuations 
are strictly correlated at all points in space-timel. This conclusion, asides 
from being highly counterintuitive, contradicts known facts about quantum 
fluctuations. A possible solution to the paradox is that correlations between 
events in different points of space cannot be written, as in eq.f l20l) , using 
joint probabilities, a well known fact in quantum mechanics (for instance, 
it is crucial in the proof of Bell's theorem [10].) There is another solution 
(which does not exclude the former), namely that Minkowski space is not 
well defined in quantized general relativity. In fact, in quantized gravity the 
metric should be quantized, meaning that the metric coefficients are oper- 
ators (see eq. (l271) below). Thus neither the distance nor the time interval 
between events are well defined. In other words, given two events of coordi- 
nates (ri, ti) and (r2, ^2) there is a quantum uncertainty about the relativistic 
interval existing between them. It is possible to state with confidence that 
two events are spatially separated if |ri — r2| >> \ti — ^2! or temporally sep- 
arated if |ri — r2| << 1^1—^2! , but it is never possible to state that they are 
light-like separated. I point out that this fact already removes the divergence 
of the two-point correlation function, shown e. g. in eq. ffTS]) . 

3 Space-time curvature due to quantum vac- 
uum fluctuations 

Working within quantized gravity the space-time structure is determined by 
the quantum state, | of the universe and the matter stress-energy tensor 
operator, T^,^ (x) , of the quantum fields at every space-time point, x. Here 
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X stands for the 4 coordinates in an appropriate reference frame, that is 

X = {Xi,X2,X3,X4} , (21) 

The study of the quantum fields in curved space-times and the gravitational 
back reaction of the fields is a difficult subject [12j. In particular the curvature 
may give rise to a modification of the vacuum stress-energy However for 
our purposes the metric is so close to Minkowskian that we may treat the 
quantum fields as if they existed in fiat space-time, although we want to 
calculate the (small) curvature induced by the vacuum fiuctuations of the 
fields. 

Our approach rests upon the existence of two quite different scales in 
the problem, namely a cosmic scale (with typical distances of megaparsecs) 
and the atomic scale of the correlations between vacuum fluctuations (which 
involves distances smaller than, say, nanometers). In the latter scale quan- 
tization is essential, but in the former we may treat everything as classical, 
as is explained in the following. For any two quantum observables, a{x) 
and b{y), at the space-time points x and y respectively, we may deflne the 
correlation 

Cab {x,y) = ($ I a (x) 6(y) I $) - ($ I a (x) I $)($ \b{y)\<!>). (22) 

Now it is plausible to assume that the correlation may be relevant at the 
atomic scale but goes to zero when the distance increases toward a macro- 
scopic scale. As a consequence, in the cosmic scale we may treat the expec- 
tations of quantum observables as classical variables, and the expectations of 
products of observables as products of the corresponding classical variables. 
For instance 

($ \a{x)b (y) I $) ~ a{x)b{y), a{x) = ($ | a (x) | $), b{y) = ($ | b{y) \ $). 

(23) 

In summary, we may ignore quantization when working with problems at 
any macroscopic scale provided we use as classical variables the expectations 
of the corresponding quantum observables. In sharp contrast, at the atomic 
scale we should work within quantized gravity. This happens in particular 
when X = y, that is 

($ I a (x) 6 (x) I $) 7^ ($ I a (x) I $) X ($ I 6 (x) | $). 

The main hypothesis of this paper is that the expectation of the stress- 
energy tensor operator of the quantum flelds at any point gives the matter 
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(baryonic or dark) stress-energy, without any additional contribution of the 
vacuum. With reference to eqs.(I3]) and (E]), this means that 

($ I f ° I $) = ($ I f; I $) ^ for /iz. ^ 00. (24) 

This suggests defining a vacuum stress-energy tensor operator as 

= %^ - I I ^)i = T,u - t;:i^'i (25) 

where / is the identity operator. The existence of vacuum fluctuations means 
that, although the expectation of TJ^^'^ is zero by definition, there are corre- 
lated vacuum fluctuations, that is 

$ |f;r (x) {y)\ $ ) ^ in general. (26) 

In order to proceed with the calculation I shall start with the quantum 
metric 

ds"^ = 'g^ydx^dx" ^ (27) 

using polar coordinates 

= t, x^ = r, = e, x^ = (f), (28) 

and I shall write the (quantum operators) coefficients of the metric in the 
form 

^00 = -1 + hoo, ?ii = 1 + hi, g22 = r'^ (l + h22^ , 

^33 = sin^ 9 (^1 + , gf.^ = h^.^ for /x 7^ z/, (29) 

(multiplication of every term times the unit operator is implicit). Here h^y 
is a (small in some sense) correction to a Minkowski metric. If we want that 
the vacuum expectation of eq.( l27|) agrees with eqs.(l3]) to ([8]) we should have, 
to order O (r^), 

V) = except (^/ioo) = {Pde + Pmat) 
where (/i^^^ stands for h^^ $^ . 
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The proof will consist of the following steps: 

1. We should define an Einstein quantum tensor operator G^,^ in terms 
of the operators 'gf^u (or what is equivalent, the operators h^i,). 

2. Assuming that in quantized gravity the counterpart of Einstein equa- 
tions reads ^ ^ 

G^u = ^ T^i,, (31) 

we should solve these (non-linear coupled partial differential) operator equa- 
tions in order to get the quantum metric coefficients in terms of inte- 
grals involving the stress-energy tensor operators TJ^ (x) and products like 

f; (x) (y) , f; (x) (y) f; {z) , etc. 

3. Finally we should calculate the expectation of the metric coefficients 
'g^u in terms of integrals involving the expectations 



r: (x) (y) $ ) , ( $ t; (x) {y) t; {z 



$ ) , etc. 



The expectation of the metric should reproduce eqs. (l30ll . 

A problem appears in the first step because there is not yet a quantum 
gravity theory specifying G^j^y in terms of which would involve a quantum 
counterpart of Riemann's theory. I will not solve the problem in general, 
but for the approximate expression of G^v containing only terms linear or 
quadratic in the (small) operators h^^, I will make the plausible assumption 
that the expression of G^^, in terms of h^,^, and their derivatives with respect 
to the coordinates, is the same as the one for the corresponding classical 
quantities with the rule that the operators should appear in symmetrical order. 
The latter assumption means that the operator corresponding to the classical 

product ab will be the quantum expression | (db + ba 



4 Quantum Einstein equation and its solu- 
tion 

In order to simplify de calculations I will introduce the approximation of 
retaining, in the expression of G^^, only terms of zeroth and first order in 
/i^i,, except for both /iqo and hu, which will be mantained up to second 
order. With these approximations our calculation simplifies substantially by 
the following reasons. Firstly it may be realized that terms of zeroth order 



13 



will not contribute to G^y{x) because to zeroth order the metric eq.( 
is Minkowskian. In addition, the terms linear in h^^ with ^ 00 and 
fiu ^ 11 (and of zeroth order in both /iqo and hu) will not contribute to the 
expectations (g^u{x)^ and (g ^u{x)G \u{y)^ when eqs.f l30|) are taken into 
account. Consequently we may ignore such terms from now on, which in 
practice is equivalent to putting /i^y = whenever jjiu ^ 00 and jjiu ^ 11. 
This amounts to replacing the metric eq.( !29|) by 

ds^ = exp (a) dr^ + (r^dO'^ + sin^ ed<f)I- exp 3) dt^, (32) 

where / is the identity operators and I have introduced the new functions 
a = log (^1 + hn^ , 3 = log (^1 - ho 



for latter convenience. Eq. fl32|) looks like the metric of a space-time with 
spherical symmetry in standard coordinates. However there are two impor- 
tant differences. Firstly the metric tensor is a quantum operator rather than 
a classical (c-number) tensor. Secondly the quantities a and (3 depend on 
the coordinates 9 and in addition to the dependence on t and r, typical of 
spherical symmetry. 

The quantized metric eq. (132|) should be used when working at the atomic 
scale, but at the cosmic scale we may use a metric obtained by the expectation 
of the former, that is 

ds^ = ($ I ds^ I $) (33) 
= ($ I exp (S) I ^)dr^ + r^de^ + sin^ Odcj)'^ - ($ | exp( /3) | ^)dt'^. 

In terms of a and /3 eqs. ( |30|) should be written as follows 

^ (Pde + Pmat) ~ ($ I exp (a) I $) - 1 ~ ($ I S + y I $), (34) 

^ (]^Pmat-PDEj ^ | exp {^) | $) - 1 ($ | 3 + y I (35) 

Our aim now is to justify these two equations as deriving from vacuum fluc- 
tuations. Thus I shall obtain the expectations of the right sides of eqs.f p^ 
and in terms of the correlations (two-point functions) of the density 
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fluctuations of tlie vacuum fields. To do that we should begin getting the 
appropriate quantum Einstein equations (involving tensor operators) and 
solving them. 

From the metric eq. ( p2|) it is straightforward to get the quantum Einstein 
equations provided we assume that they are similar to the classical counter- 
parts, as explained above. Two of them do not contain time derivatives and 
they are the only ones to be studied here, that is 



8nGp = G° = ^- 


1 doL 


a da 
r Or 
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2r2 




1 d^a 


1 d'^a 






1 


2r2 de^ 


2r2s2 d(i? 


¥2 ( 


'39 ) 





(36) 



-8.Gp = G} = 4-4--^ + -^-Acot^^ (37) 



r 



2^2 r dr r dr 2r2 06 



1 d'(3 1 d'/3 1 fdpy 1 {dp 



2r^ 86^ 2r2s2 50^ Ar^ \d9 J ir'^s'^ \d(p 

Here I have removed the carets of the operators for notational simplicity 
which will be also made from now on. But I remember that both a, /3 and 
their derivatives are quantum operators and that whenever we have a product 
of two of them it is understood that it means symmetrically orderer product. 
For instance 

ad/3 „ 1 fad'f3 d^a 

— — actually means - — — h 

r or 2 \r or or r 

After some algebra eqs.( l36l) and ( 1371) may be rewritten, in more compact 
form, 

9^ da da 1 ^ 1 ^ 1 * / 9x 

8Txr^Gp = a + r— - ra- Aa + -aAa - -/^ia^), (38) 

2 or or 2 4 8 

- ^-r^Gp = « - Y - + ™f - + i/3A/3 - ^A , (39) 

where A is the angular part of the Laplacian operator, that is 

I d { . 1 
A = ^-TTT^ sm6'— + 



smOde V dej sm'ed(f) 
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I shall start solving the nonlinear partial differential eq.f l38p , which con- 
tains a single unknown, namely the operator a {r, 6, 0) . Actually a may also 
depend on time, but as eqs.fl36l) and fl371) do not contain time derivatives 
the time t appears as a parameter, rather than one of the variables of the 
partial differential equations, and I shall not write it explicitly. In order to 
solve eq.f lSS]) I will approximate the solution by a perturbation expansion 
in powers of the Newton constant G, and work to order O(G^), writing the 
(operator) metric parameter a in the form 

a = Gao + G^ai, (40) 

If I put this in eq. flSSl) the terms of first order in G give the linear (in the 
unknown ao) equation 

n dan 1 , , , , 

87rrV = «o + r^--A«o, 41 
or 2 

whilst the terms of second order give an equation also linear (in ai), namely 

ai + r^ - -Aai = ^ + ra^—^ - -aoAao + -A (al) . 42 
or 2 2 or 4 8 

The solution of eq.( l4T]) , with the condition that ao = at r = 0, may be 

written in simplified notation as 

ao = 87rA(rV), (43) 

with the meaning 

ao (r) = 87r / dh^A (r, r^) [rjp (n)] , r = {r, 9, 0} , (44) 

^0 



where A is a kernel to be specified later on. Eq.( l42l) , with the condition that 
ai = at r = 0, may be solved similarly leading to 



ai = A 



an dao 1 A 1 A / 2\ 

-TT + TOo^ TttoAao + - A (aoj 

2 dr 4 8 



(45) 



where ao is given by eq.(H3l) . Thus the solution of eq.( !38|) may be written, in 
simplified notation (I shall use units G = 1 from now on, although writing 
explicitly Newton's constant sometimes for the sake of clarity). 



a = 87rA(rV)+647rM 



^ ( 1 + r|- + Ia) [A{r'p)r - \[A{r'p)]A[A{r'p)] 

(46 
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We are interested in the expectations eqs.(lMI) and we get 

Q? O? O? 

($ I a + y I $) = ($ I a + y I ^)^at + ($ | « + y | $)^ac, (47) 



($ I a + y I ^)„,at = 87rA(rV„aJ + 327r2A 



4 + 7A ) [A{r'Prnat)f 



dr 4 

IQ-K^A [[A(rV„.at)]A[^(rV™at)]] + 327r2[A(rV^aJ]', (48) 



(<f I « + y |<|.)_^32vr2A(l+r|- + iA)($|KrV_)]^|$), 

-167r2($ I [A(rV.J]A[A(rV.J] |$)+ 327r2(<|.|[A(rV.J]'|$), (49) 

where I have taken into account eqs.f l24l) and (!25l) . The proof is not difficult 
taking into account that we work to ffist order in ai and to second order in 
Uq. Let us consider for instance the term 

($ I «^ I $) ~ ($ I «o I $) 

= Ut:^ j dhirlA{r,r^) j rf^rar^A (r, r^) ($ |p (n) p (rs)! $). 

Taken eqs.f E^ and f l2^ into account the two-point correlation becomes 
($ |p (ri) p (ra)! $) = ($ | [p„,t (ri) + p,„, (ri)] [p^„, (rs) + p,,, (ra)] | $) 

A similar analysis may be made for the other terms. The result is that the 
expectation ($ | a + a^/2 | $) is the sum of two expressions, one containing 
only the matter density, PmatJ '^''^^ other one the vacuum density, p^^^, 
i. e. there are no cross terms. It may he also realized that we should solve 
eg. fl38l) at least up to terms in in order to get the leading term due to 
the vacuum fluctuations. This is because the expectation of the solution to 
order G, eq.flHl) , gives no contribution due to the vanishing of {vac |p| vac) , 
see eq.(I21]) • 

The terms with p^^^ give the contribution of matter density to the metric 
coefficient a. In particular, if we model the matter density of the universe by 
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a constant it is not difficult to check, taking eq.f l54p into account (see below), 
that the matter term gives 

$ a + -a^ ~ + M = —p^atr\ 50 

2 r r-^ 6 

which agrees with the second order expansion of exp a in the well known 
Schwarzschild solution 

exp a = I 1 I 

5 Contribution of the vacuum fluctuations to 
the metric 

In the following I shall calculate the different terms involved in eqs. (l49l) . 
I start solving eq.lHTl) by writing p (r, 6*, 0) and ao{r,9,(j)) as expansions in 
terms of spherical harmonics, that is 

Im ^ 

ao = '^aimir)Yimi9,(f)) ^ atmir) = j ao{r, 9, (f))Yi'^ (6, (p) dn. (51) 

Im 

1 get from eq.(HTl) 



^T^Plmr^ 



aim + r—j—, (52) 



whose solution with the initial condition aim (0) = is 

aim (r) = Svrr f {y /rf+'^'^'^" p,^ (y) dy, (53) 
Taking eqs.( !5T]) into account I get (see eq.( H4|) ) 

^(r, n) ^ r-' J2 Yirn {e, 0) in/rf'^'y^Y;^ (6,, <P,) . (54) 



Im 
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I am interested in the expectations involving p^^^ defined in eq. (jlSD , where 
I shall write (.) for ($ |.| $) for notational simplicity. I begin with 

X dr2 j d^2riA{v,V2)C{s), (55) 

where the two-point correlation function of the density, C(s), was given in 
eq.(|T2|) with 



s = |ri — = y fi + ^2 ~ 2rir2M, u = cos ^12. (56) 

A consequence of the assumption that the correlation C depends only on the 
distance |ri — r2| , rather than on ri and r2 separately, is that the expectation 
eq.( l55|) will not depend on the angular variables {6,(f)) and we may average 
over these variables. Hence it follows that I' = l,m' = m and we get 



{{A[r^p]f) = J dnjdn.j dr2jdn2 

J2(^ir2/rY'+'y'^'Yi^ {01, 0i) {02, 02) C{s). (57) 



Im 

This leads to 

„2 



„2 i\2\ ^ 



f dri f dr2 f du V {21 + 1) {nr2/r^ f+'^'+^^'^ Pi {u) C{s 
Jo Jo J-1 I 

(58) 

where I have taken into account the following property of spherical harmonic 
functions 

Yin. (^1, 0i) Y;^ (^2, 02) = (cos 9^2) , (59) 

m 

6*12 being the angle between the directions (6*1,0]^) and (^^2,02) and Piiu) a 
Legendre polynomial. Now I introduce the new variables {y, z} defined by 

ri = y + z/2, r2 = y - z/2, (60) 
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substitute an s integration for the u integration and change the order of the 
integrals. Thus eq.( 155]l becomes after some algebra 



(a 



-1 roo ps 

'■) ^ {(A[r'p])')^—j^ C{s)sdsY,{'^l + ^) j^dzk{r,z) (61) 



r-z/2 



dy 



z/2 



y - z- 



74 Y 



Pi I 



2y2 - ^2/2 J ■ 



A similar method allows calculating the term {^\{Ap^^JA{Ap^^J\^). 
Putting all relevant terms of eq.f l49p together we obtain 



(a) = SttA 



C(s)st^s^(2/ + 1) J dz i^l 



(62) 

It may be realized that the main contribution to (a) comes from high values 
of I. Thus we might neglect 1 <^ |Z (Z + 1), which also shows that the term 
is negligible in comparison with the term — ^aoAao in eq.f HS]) . Similarly 
we may neglect 



< (a) 



(63) 



in eq. ( 149|) . Proceeding with the calculation of eq. (162|) , I get after some alge- 
bra 



l + r^ + lni + l)]liir,z) 



r — z 



l+/(/+l)/2 



P, 1 



2r2 



! + -/(/ + !) 



h {r, z) , 



(64) 



where I have neglected /2 <^ 2r'^ in the argument of the Legendre polyno- 
mial. In the quantity Ii (r, z) we may approximate 2y'^ — z'^ /2 by 2r'^ because 
only values of y close to r contribute substantially to the y-integral. This 
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leads to 
// ~ P/ ( 1 - 



Plil 



2r2 
s'^ — 



z/2 \ 



z 



z 



exp -— [2 + /(/ + 1)] 



[2 + /(/ + !)] 



zl2 



r[3 + /(/ + ir'Pi 1 



2r2 



X exp -— [2 + l{l + 1)] 



/r,\ 2+Z{«+l) 

r — zj 1\ ' 



.2 ^2- 



r [3 + / (/ + 1)]-^ ( 1 - '^-^^ exp (-|: [2 + /(/ + 1)] 



(65) 



where I have neglected 



|:^[2 + /(/ + i)]«|:[2 + /(/ + i)] 



in the exponential because z < s ^ r. Hence eq.flMl) becomes 

4 + /(/ + !) 



(a) ~ 47rA 



/■oo 

/ C(s)srfs^(2/ + 1) 



3 + /(/ + 1) 



X / d^exp(^-— [2 + /(/ + !)] jP; (1 



2r2 



(66) 



We may substitute a Bessel function for the Legendre polynomial, which is 
a good approximation for large / and argument close to unity, as is shown by 
the well known limit (see Gradshteyn et al.[T3], n" 8.722-2) 



lim Pi (cos Y ) = Jo (x) . 



After that I will substitute an integral for the sum in / and I get 



/'CO /"S /'OO 

47rA r / C{s)sds dz Idl exp [-zP /2r] Jq 
. Jo Jo Jo 



AirA 



r dz 

I C{s)sds / — exp 
Jo ^ 



2zr 



(67) 
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where the / integral has been taken from the hterature (see Gradshteyn[ 
n" 6.631-4) . The z integral may be performed in terms of the new variable 
t = s^/(2rz), which leads to 

o2 



dz 
— exp 

z 



z^ — 



2rz 



s/2r 
oo 

s/2r 



dt 




T 


exp 


dt 




T 


exp 



2rH 



t 



s/2r 



dt . , 

— exp e exp [—t\ 



exp [-t] = log(2r/s) - 0.557, 



(68) 



where I have put expe ~ 1 because e = j {2r'^t) <^ 1 for t > s/2r. Now, 
taking into account the action of the kernel A, eq.( l54l) , applied to a function 
/ (r) not depending on the polar angles {9, 0}, we get finally 

„2 



911 I $)mc ^{<^\a + 



a 



(69) 



C{s)sds (log(2r/s) - 0.557) ~ GOOr'^ / C{s)sds 



where I have taken into account eg. (1631) and I have estimated r/s ^ 10'^°. If 
we include the contribution of matter and put explicitly Newton's constant 
we get 

SnG 



a 



($ I a H I $) 



-p_,r^ + 600GV / Cis)sds. 



(70) 



Now I shall solve the second component of the Einstein equation, eq.( l39|l , 
which involves both coefficients a and /3 of the metric. I will search for a 
solution of the form 



(3{r,9,4>) = -a{r,9,(f))+-f{ 
which taking eq.f l38|) into account leads to 



+ O (r^) 



Svrr [p + p) = 2r 7 — 2r 07 + -r A7 — Aa + -aAa 

-^{a - jr')A{a - ^r^) - \A{a^) + ^A[(« - ^ir^f] 

= 2r^j -Aa + (r^) . 

where in the second equality I have taken into account that a is of order 
O (r^) (see eqs.fjSTll and (j53ll). This leads to the solution 

(3 = -a + 4nr^ {p + p) + ^ Aa + O (r*^) . 
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Hence eqs.f l34l) and fl35p lead to the following equalities, where terms of order 
O (r^) are neglected 



^2 



mat 



/3 



mat ■ 



The third equality takes into account that ($ | S | $) does not depend on the 
angles ^,0, so that ($ | AS | $) = A($ | S | $) = 0, and the last equality 
the fact that {3 is of order O (r^) . Hence, taking eq.( !70|) into account, we get 



I /3 + Y I ^) ^ -^Pmatr' - Cis)sds. (71) 

6 Conclusions 

The main result of the paper is that eqs.f lTU]) and fITT]) provide the expectation 
of the quantized metric (see eq.( !33|) ). in terms of the two-point correlation 
of the vacuum fluctuations, C(s), that is 

sin^ ed(P^ - (^^p^y - GOOGV^ C{s)sd^ dt\ (72) 

This is identical to the standard FRW metric, eq.(IT]) , via the approximate 
"free falling" metric eq.dH]) , provided that we identify 

Pjje ^ WOG / C{s)sds, (73) 
Jo 

which shows that vacuum fluctuations give rise to a curvature of space-time 
similar to what would be produced by a "dark energy" density (plus cold 
matter). However we cannot fix the value of pj^^ as far as we do not know 
the two-point correlation function of vacuum density fluctuations. Crucial 
for the result is the fact that Einstein's equation involves a non-linear re- 
lation between the stress-energy tensor and the metric tensor. In fact, if 
the relation was linear, then the vanishing of the vacuum expectation of the 
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quantum matter density operator would imply vanishing of curvature, that 
is Minkowski space, in the absence of matter. 

Our results suggest that the observed accelerated expansion of the uni- 
verse might be explained as due to the quantum vacuum fluctuations, without 
the need of any "dark energy". If this is the case dark energy, p^^;, appears 
as a parameter which mimics the effect of the vacuum fluctuations. The de- 
rived relation between the two-point correlation of the fluctuations and the 
value of the dark energy parameter would allow calculating the latter if the 
said two-point correlation was known. Arguments are given which suggest 
that such calculation might give results in agreement with observations. 



7 Appendix. Two-point correlation of vac- 
uum fluctations of the radiation field. 



Let us start calculating the two-point correlation 
A{ri,ti.r2,t2) = (vac 



:E(ri,ti)':)(:E(r2,t2)' 



vac 



(74) 



Taking eqs.f ll7p into account it is not difficult to prove that contributions 
to A, eq.( !74|) , will derive only from terms with two annihilation opera- 
tors coming from (: E(ri,ti)^ :) and two creation operators coming from 
(: E (r2,t2)^ 0- Thus the contributing terms will derive from 

— V ww'akeCtk'e'^ (k) ■ s' (k') exp [i (k' + k) .ri — iuti] 



ke k'e' 



fl 



X ■ 



V 



/"a+,„a+,,,„e"(k")-s'"(k'") 



k"e" k"'6"' 

X exp [-i (k" + k'") .rs + ioot2] . 
Taking into account that 



[vac \a\ieO']i'e'a^"e"'^t"'e"' \ '^^^) = <^kk"'^k'k"' + 5kk""^k'k", 



we get 
A 



2f^ 



2fi 



ke k'e' 



k k' 



(k • k'y 

k2k'2 



exp [i (k' + k) .r-i {u + u') t] (75) 
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where r = ri— r2,t = ti — t2 and I have performed the sum in polarizations 
in the second equahty. I shall use polar angles taking the direction of the 
vector r as polar axis, so that 

k =u (sin 9 cos 0, sin 6 sin 0, cos 9) , k =uj' (sin 9' cos 0', sin 9' sin 0', cos 9') . 

After the standard replacements 



r2n 

d(j), 



eq.f lTS]) gives 

A = —-^ tu^dco w'W / d{cos9') d<p' d{cos9) di 

rJ t/Q — 1 •-'0 ^ — 1 •-'0 

X exp [iu' (r cos 9' — t)] exp [iu' (r cos 9' — t)] 

X 1 + (sin 9 cos sin 9' cos 0' + sin ^ sin sin ^' sin 0' + cos 9 cos 0')^ 
Integrating the angles and 0' we get 

^2 rcxD i>co i>l i>l 

A = J u^dco J u'^du' J du J du[3 + 3u'^u'^ - - u'^\ 

X exp [iw [ru — t)] exp [iw' [ru — t)] , 

where r = |r| , m = cos9,u' = cos 9'. The u and w' integrals are trivial and I 
obtain 

^2 roo r-oo 

A = J / oo^duj / w'^rfw' exp \-i {uj + uj') t] 

16vr^ Jo Jo 

X [3/ (x) J (x') + 3 J (x) J (x') -I{x)J (x') -J{x)I (x')] , (76) 



where 



/ (x) = / exp (— ixw) = — ^- — —, X = ur, 

J-l X 

(i^ 2 sin X 4 cos x 4 sin x 

W = -T^^W = — + 

dx^ X x^ x"* 
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Now 

r . • ^ w ^ '2 r 2, . • / ^ . 2 (3t^ + r^) 

/ u auexni—tutjllx) = —r x aujexp(—txt/r)smx = 5-, 

Jo V ^ V ^ ^4 P I / J _ ^2^3 

co^ dco exp (—icot) J (x) = — co^ du exp {—ixt / r) -— i 

2 {t^ + 3r^) 

The correlation between the magnetic energies is the same as eq.f TMl) .The 
correlation between electric and magnetic parts may be derived without dif- 
ficulty. It is 

Of) 2 

B = -^^^a;a;'[z(krX£(k)) ■6'(k')]'exp[?(k' + k).r-z(a; + u;')^], 

ke k'e' 

k.r 

kr = — =u; (0, 0, cos6') , 



whence 



B = -— coco' (1 + cos^ 9) exp [i (k' + k) .r-i {co + co') t] 



ke k'e' 

^2 roo 



in 



^0 



^ / (^^du / u'^du' (77) 



X [/ (x) / (a;') + J{x)J {x')] exp [-i (u + co') t] . (78) 
Taking into account eqs.(l76l) to ( l79l) we finally obtain 

C{r,t) = 2A + B= (79) 
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